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SUMMARY 



Understanding BRCAI mutant cancers is hampered by difficulties in obtaining primary cells from patients. We therefore generated and 
characterized 24 induced pluripotent stem cell (iPSC) lines from fibroblasts of eight individuals from a BRCAI 5382insC mutant family. 
All BRCAI 5382insC heterozygous fibroblasts, iPSCs, and teratomas maintained equivalent expression of both wild-type and mutant 
BRCAI transcripts. Although no difference in differentiation capacity was observed between BRCAI wild-type and mutant iPSCs, there 
was elevated protein kinase C-theta (PKC-theta) in BRCAI mutant iPSCs. Cancer cell lines with BRCAI mutations and hormone-receptor- 
negative breast cancers also displayed elevated PKC-theta. Genome sequencing of the 24 iPSC lines showed a similar frequency of reprog- 
ramming-associated de novo mutations in BRCAI mutant and wild-type iPSCs. These data indicate that iPSC lines can be derived from 
BRCAI mutant fibroblasts to study the effects of the mutation on gene expression and genome stability. 



INTRODUCTION 

BRCAl/2 mutations are autosomal-dominant mutations 
that dramatically increase the risk for developing breast 
and ovarian cancer, and to a lesser extent other malig- 
nancies, such as melanoma, pancreatic cancer, and pros- 
tate cancer (Futreal et al., 1994; Lancaster et al., 1996; 
Miki et al., 1994). The BRCAI and BRCA2 tumor-suppres- 
sor proteins play roles in transcriptional regulation and 
DNA repair (Turner et al., 2004). BRCAI also participates 
in cell-cycle regulation (McPherson et al., 2004), polyade- 
nylation of messenger RNA (mRNA) (Kleiman et al., 
2005), and ubiquitinylation (Baer and Ludwig, 2002). The 
cancers that arise in patients with inherited BRCAI muta- 
tions appear to be more aggressive than those in patients 
with BRCA2 mutations or sporadic breast tumors. This 
aggressiveness is thought to result at least in part from 
the fact that Bi^C^ J -deficient tumors, but not BRCA2-deh- 
cient tumors, are usually estrogen receptor (ER) and proges- 
terone receptor (PR) negative (Turner et al., 2004) and 
therefore unresponsive to hormonal therapies. These 
observations suggest that the shared DNA repair functions 
of BRCAI and BRCA2 do not solely account for differences 
in cancer incidence among patients with mutations in 
these genes. Evidence also indicates that cancer risk differs 
depending on the location of the mutation in the BRCA 
sequences (Risch et al., 2006). Examination of different 
BRCA 1/2 mutant cells that are prone to tumor formation 



will identify new BRCAI and BRCA2 tumor-suppressor 
functions. 

One BRCA2 {BRCA2 6174 delT) and two BRCAI (BRCAI 
5382insC and BRCAI 185 delAG) founder mutations 
account for 90% of inherited BRCA 1/2 mutations in 
patients of Ashkenazi Jewish descent (Petrucelli et al., 
2010). Because 0.1% of the Ashkenazi population carries 
the BRCAI 5382insC mutation, we estimate that this muta- 
tion is present in approximately 10,000 individuals in the 
United States alone. A deeper understanding of how each 
of the Ashkenazi BRCA mutations leads to breast cancer is 
important because it will improve our understanding of 
phenotypic variabilities due to different mutations in the 
same gene. 

Although mice genetically engineered to have either 
Brcal or Brca2 deficiency have been informative models 
for studies of breast cancer development (Drost et al., 
2011; Drost and Jonkers, 2009; Evers and Jonkers, 2006; 
Shakya et al., 201 1), differences in underlying biology exist 
between humans and mice. Thus, human models are 
necessary to complement animal models. It has been re- 
ported that mutation of a single BRCAI allele leads to 
genomic instability in human cells, a phenomenon not 
observed in mice (Konishi et al., 2011). Generation of 
induced pluripotent stem cells (iPSCs) (Takahashi et al., 
2007) from patients carrying BRCAl/2 mutations may pro- 
vide a window into the cellular phenotypic differences 
driving their increased cancer risk. BRCAl/2 mutant iPSC 
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Figure 1. Pedigree and iPSC Generation 
and Characterization Flowchart 

(A) The four-generation pedigree of the 
BRCAl 5382insC family. Patient identifiers 
are Located beLow each symbol. Patient ages 
at the time of sample acquisition or death 
are Located to the upper Left of each symboL. 
Shaded symbols identify patients carrying 
the BRCAl 5382insC mutation. A diagonal 
line indicates deceased individuals. Aster- 
isks identify patients whose skin samples 
were used for iPSC derivation. 

(B) Flowchart for iPSC generation using 
nonintegrative mRNA and miRNA-based 
reprogramming, gene expression by micro- 
array, qPCR, western blots, and whole- 
genome sequencing. 



lines may also prove useful for screening therapeutic com- 
pounds for the prevention and treatment of genetically 
predisposed cancer patients. 

To probe the cellular phenotypes of patients with BRCAl- 
associated cancers, we reprogrammed primary dermal fi- 
broblasts from a family carrying the Ashkenazi BRCAl 
5382insC mutation into 24 iPSCs. Here, we report on our 
characterization of these cell lines using in vitro and in vivo 
cellular differentiation assays as well as gene-expression 
and whole-genome sequencing (WGS) analyses. 



RESULTS 

Generation of Wild-Type and BRCAl 5382insC 
Heterozygous iPSCs 

Prior to this work, it was not known whether iPSCs could 
be generated from patients with BRCAl mutations, due 



to the reported genome instability of cells heterozygous 
for these mutations (Konishi et al., 2011). We recruited a 
large family with BRCAl mutant cells to donate samples 
for fibroblast reprogramming (Figure lA; Table 1). We ob- 
tained skin biopsies from family members who spanned 
three generations and ranged in age from 20 to 90 years. 
Four family members had no mutation at the BRCAl locus, 
and four were heterozygous for the BRCAl 5382insC muta- 
tion (Figure lA). A prior reference to this family (Ross, 
2012) noted that the index breast cancer patient (111.5) 
died at age 38. This patient was heterozygous for the fam- 
ily's BRCAl 5382insC mutation, as her 22-year-old 
daughter (IV.3) carries the mutation (Figure lA; Ross, 
2012). 

To avoid the complications of insertional mutagenesis, 
which often occurs with retroviral expression of reprog- 
ramming factors in cells from an already cancer-prone 
BRCAl family, were generated iPSC lines from fibroblasts 
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Table 1. Clinical Characteristics and BRCAI Status of the Study Subjects 


Patient No. 


Age° and Gender 


Cancer History (Age at Diagnosis) 


BRCAl'' 


Other Health Issues 


II. 1 


80 M 


none 


WT 


none 


II.2 


87 M 


nanrrpatir ranrpr ^85^ rninn ranrer ^70^ 


WT 


none 


II. 3 


65 F 


nvarian ranrpr ^64^ 

UVQIIQII V-Oll^d \\J^I 


WT 


none 


II.4 


60 M 


aHrpnal ranrpr 


WT 


none 


II.5 


83 M 


pheochromocytoma (50), lung adenocarcinoma 

^RS^ mplsnnma hrn nrhna 1 \/pnl 3 r 1 1 1 nn ^Rfl^ 


WT 


none 


II.6 


90 F 


pnrinmptrial ranrpr (^()\ hrpa«;t ranrpr ^7'^^ 


mutant 


mnnnrlnnal nammnnathu nf iinHptprminpri 

significance 


III.l 


50 M 


none 


WT 


dysplastic nevi, chondroma (elbow), colon polyps 


III.3 


52 M 


testicular cancer (28), prostate cancer (55) 


mutant 


none 


III.4 


45 F 


none 


WT 


none 


III.7 


48 F 


melanoma (38) 


mutant 


dysplastic nevi, endometrial atypia, focal nodular 
hyperplasia (liver), prophylactic surgeries' 


IV. 1 


19 F 


none 


WT 


none 


IV.3 


22 F 


none 


mutant 


none 


''Age at the time of sample acquisition or at the time of death. 

''WT, wild-type for both BRCAI alleles; mutant, BRCAI 5382insC heterozygous. 

"^Bilateral mastectomies, hysterectomy, and oophorectomies. 



using nonintegrating mRNA-based transfections to ex- 
press the reprogramming factors (0CT4, S0X2, KLF4, 
LIN28, and c-MYC) (Warren et al., 2010). Primary fibro- 
blasts were isolated from skin biopsies (II. 1, II. 6, Ill.l, 
III.3, III.4, III.7, IV.l, and IV.3; Figure IB), and at passage 
5 or 6, these cells were transfected with reprogramming 
factors and GFP (Warren et al., 2010). All eight fibroblast 
cultures exhibited similar transfection efficiencies, as 
measured by GFP protein expression during the course of 
daily transfections (data not shown). The majority of 
patient fibroblasts were reprogrammed within two at- 
tempts. Two samples, II. 1 and lll.l, both of which are 
wild-type at the BRCAI locus, failed to be reprogrammed 
by mRNA alone after three attempts. They were then sub- 
jected to either one of two cocktails of microRNAs 
(miRNAs) previously used for reprogramming by Anokye- 
Danso et al. (2011) (mixture A) and Miyoshi et al. (2011) 
(mixture B). Cotransfection on days 1 and 4 with either 
miRNA mixture A or mixture B with the five-factor 
mRNA transfections led to the generation of iPSC-like col- 
onies. We did not observe a significant change in iPSC 
colony frequency (colonies per number of cells trans- 
fected) in the presence of the BRCAI 5382insC mutation 
(Table SI available online). On average, three indepen- 
dently derived clonal iPSC lines were selected and 
expanded from each patient's fibroblasts (Table 2). 



Differentiation of Wild-Type and BRCAI S382insC 
iPSCs 

The 24 iPSC lines were characterized for pluripotency. As 
illustrated in the representative examples in Figure 2A, all 
isolated iPSC lines displayed an embryonic stem cell-like 
morphology and expressed pluripotency-associated 
antigens including TRA-1-60 (proteoglycan), TRA-1-81 
(proteoglycan), SSEA3 (carbohydrate), SSEA4 (glycolipid), 
and NANOG (Figures 2A and SI; data not shown). No 
significant differences in morphology or other growth 
patterns between the BRCAI wild-type and mutant iPSCs 
were noted. To further evaluate the pluripotency of the 
iPSCs and to determine whether functional differences in 
differentiation capacity existed between the wild-type 
and mutant cell lines, all iPSC lines were induced in culture 
to differentiate into cells representative of all three germ 
layers. Analysis of lineage-specific protein expression 
by immunofluorescence indicated that the iPSC lines 
readily differentiated into cells of the ectoderm (assessed 
by PAX6 and pill-tubulin expression), mesoderm 
(BRACHYURY and GATA4 expression), and endoderm 
(FOXA2 and SOX17 expression; Figures 2B and S2). No dif- 
ferences in differentiation potential between wild-type and 
BRCAI 5382insC mutant iPSCs were observed. 

To observe in vivo differentiation capacities of BRCAI 
wild-type and mutant iPSCs, each cell line was injected 
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Table 2. Summary 


of iPSC Derivation from Fibroblasts 






Fibroblast Samples' 


Reprogramming Method'' 


iPSC Clones 


Attempts to Reprogram BRCAl Status*" 


II.l 


mRNA + mlRNA 


1PS.II. 1.1, IPS. II. 1.2, IPS. II. 1.3, IPS. II. 1.4 


4 


WT 


II.6 


mRNA 


1PS.II. 6.1, IPS. II. 6. 2, 1PS.II. 6. 3 


1 


mutant 


III.l 


mRNA + mlRNA 


1PS.III.1.1, 1PS.III.1.2 


4 


WT 


III.3 


mRNA 


1PS.III.3.1, 1PS.III.3.2, 1PS.III.3.3 


2 


mutant 


III.4 


mRNA 


1PS.III.4.1, 1PS.III.4.2 


1 


WT 


III.7 


mRNA 


iPS.III.7.1, 1PS.III.7.2, 1PS.III.7.3, 1PS.III.7.4 


2 


mutant 


IV.l 


mRNA 


1PS.IV.1.1, iPS.IV.1.2, iPS.IV.1.3 


2 


WT 


IV.3 


mRNA 


1PS.IV.3.1, 1PS.IV.3.2, 1PS.IV.3.3 


1 


mutant 


See also Table SI. 

'Sk1n samples were harvested from patients after they gave Informed consent. 

''mRNA transfection (Warren et al., 2010) and miRNA choices and expression were as previously described (Anokye-Danso et al. 
"WT, wild-type at both BRCAl alleles; mutant, BRCAl 5382insC heterozygous. 


, 2011; Miyoshi etaL, 2011). 



into the kidney capsule of nonobese diabetic (NOD)/severe 
combined immunodeficiency (SCID) mice and observed 
for teratoma formation and composition. Once tumors 
formed, the histopathology of the tumors was evaluated 
for mesoderm-, ectoderm-, and endoderm-derived tissues. 
The ability to form benign teratomas containing tissues 
from all three germ layers is an inherent characteristic of 
iPSCs. Figure 2C displays histologic sections from represen- 
tative BRCAl -mutant- and wild-type iPSC-derived tera- 
tomas, and demonstrates the presence of tissues from all 
three germ layers. The ectoderm is represented by the neu- 
ral epithelium, the mesoderm is represented by muscle, 
and the endoderm is represented by gut epithelium. All 
iPSC lines formed teratomas in vivo without obvious differ- 
ences in their differentiation patterns (Figure S3). Invasive 
cancers were not observed in any of the BRCAl wild-type or 
mutant teratomas. 

BRCAl mRNA Expression Levels in Fibroblasts, iPSCs, 
and Teratomas 

Because the BRCAl mutant fibroblasts reprogrammed so 
readily, we hypothesized that rare fibroblasts with loss of 
heterozygosity (LOH) at the BRCAl locus or loss of 
BRCAl expression might have been present in the skin 
biopsies, and that these cells were under positive selection 
during reprogramming. The process of reprogramming has 
been reported to be similar to tumorigenesis (Riggs et al., 
2013), which is supported by the fact that p53 deficiency 
enhances reprogramming (Hong et al., 2009). To compare 
total BRCAl expression levels in fibroblasts and iPSCs, we 
measured the absolute levels of BRCAl transcripts using 
quantitative PCR (qPCR). Wild-type and heterozygous 
mutant fibroblasts expressed, on average, similar levels of 
total BRCAl with 70 copies/ng total RNA, and a 2-fold up- 



regulation of total BRCAl transcript was observed in all 
derivative iPSCs (n = 24) upon reprogramming (Figure 3A). 
This differential expression of the BRCAl gene was 
confirmed at the protein level (Figure 3 A, inset). Because 
BRCAl is thought to be involved in differentiation of 
luminal mammary progenitors into luminal epithelial cells 
(Lim et al., 2009), it was somewhat of a surprise that BRCAl 
levels increased upon dedifferentiation; however, this in- 
crease may be required to maintain genomic stability in 
iPSCs. In fact, increased expression of DNA damage repair 
genes in two iPSC lines as compared with two differenti- 
ated fibroblast lines have been observed previously (Mom- 
cilovic et al., 2010). 

To explore the possibility of BRCAl LOH or a compensa- 
tory alteration in the expression levels of either the BRCAl 
wild-type or mutant alleles in the mutant iPSCs during 
reprogramming, we quantified Bi^C^i wild-type or mutant 
allele-specific expression in the patient fibroblasts, iPSCs, 
and teratomas. To quantify the absolute levels of BRCAl 
wild-type and mutant transcripts, we designed mis- 
matched primers that specifically amplified either the 
wild- type or mutant mRNA (Figure SI). Wild-type and 
mutant BRCAl transcripts were present at equivalent levels 
in the heterozygous fibroblasts (Figure 3B). Thus, mutant 
cells (II.6, III.3, III. 7, and IV.3) exhibited a 50% decrease 
in the absolute amount of wild-type BRCAl transcript 
compared with fibroblasts with two wild-type BRCAl 
alleles (II.l, III.l, IIL4, and IV.l). All of the mutant iPSC 
lines (Figure 3C) and teratomas (Figure 3D) also expressed 
equivalent levels of both wild-type and mutant transcripts, 
indicating that neither the wild-type nor the mutant alleles 
were lost during the reprogramming process or during dif- 
ferentiation into teratomas, and that BRCAl haploinsuffi- 
ciency persisted. Hence, there was no evidence for 
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Figure 2. Generation and Analysis of iPSC 
Lines from Wild-Type and BRCAI/ 
5382insC Heterozygous Fibroblasts 

(A) Representative iPSC colonies derived 
from BRCAI wiLd-type (III.4.2) and het- 
erozygous (II. 6. 1) mutant fibroblasts. All 
iPSC lines stained positive for the pluripo- 
tency markers TRA-1-60, SSEA3, and 
NANOG. The nucleus was stained with DAPI. 
The scale bar represents 100 i^m. 

(B) In-vitro-directed differentiation and 
analyses of lineage-specific genes in repre- 
sentative wild-type and mutant iPSC lines 
from (A). The ectoderm was identified with 
PAX6 and P-tubulin, mesoderm was identi- 
fied with BRACHYURY and GATA4, and 
endoderm was identified with F0XA2 and 
S0X17. The nucleus was stained with DAPI. 
The scale bar represents 50 |im. 

(C) Representative histological analysis 
showing in vivo differentiation of wiLd-type 
(WT, top) and mutant (bottom) iPSC lines 
into teratomas with tissue components that 
represent all three germ layers: ectoderm 
(neural epithelium; patients III. 4. 2 and 

III. 3. 2), mesoderm (muscle; patients IV. 1.3 
and II. 6. 2), and endoderm (gut; patients 

IV. 1.3 and III. 7.1). Sections were stained 
with hematoxylin and eosin. The scale bar 
represents 50 (im. 

See also Figures S1-S3 and Table SI. 



silencing or loss of the wild-type BRCAI allele during re- 
programming. This was a surprise because LOH of tumor 
suppressors, such as p53, promotes reprogramming 
(Hong et al., 2009). 

Transcriptional Profiles of Mutant and Wild-Type 
iPSCs 

Since the ability to differentiate in vitro and in vivo was 
similar between BRCAI wild-type and mutant iPSCs, we 
evaluated the mutant and wild-type iPSCs for molecular 
differences. We used genome-wide microarray expression 
analysis to determine whether the BRCAI mutant iPSCs 
had a unique gene-expression signature compared with 
the wild-type iPSC lines. Expression differences (> 2-fold) 
between the mutant and wild-type lines were evaluated 
further only if those expression changes were observed in 
all mutant lines. Only one gene, PRKCQ, the gene that en- 



codes protein kinase C-theta (PKC-theta), was dysregulated 
in all of the 13 mutant iPSC lines compared with the 11 
wild-type lines (Figure 4A). No differential gene expression 
between wild-type and mutant iPSC lines was observed for 
other members of the PRKC gene family (Figure 4B). The 
specific upregulation of PKC-theta in BRCAI mutant iPSCs 
was readily confirmed by qPCR (Figure 4C) and western 
blot analysis (Figure 4D). Analysis for differential expres- 
sion of PKC-theta in BRCAI wild-type and mutant fibro- 
blasts was not possible because PKC-theta expression was 
undetectable in the fibroblasts by qPCR (Figure 4C). 
Notably, gene-expression analysis of the eight dermal fibro- 
blast cultures did not uncover any specific gene-expression 
changes between BRCAI mutant and wild-type fibroblasts. 
As might be expected for two different cell types, compari- 
son of gene expression between fibroblasts and iPSC lines 
uncovered thousands of differentially expressed genes. 
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with approximately 5,000 upregulated and 5,000 downre- 
gulated genes (data not shown). 

PKC-theta Protein Expression in Breast Cancer Tissue 

These data document the association of elevated PKC-theta 
levels with a BRCAl mutation. Because the overexpression 
was readily detected in dedifferentiated iPSCs, we sought 
to examine PKC-theta overexpression in breast tumors, 
another dedifferentiation condition. PKC-theta is a mem- 
ber of the tumor-promoting PKC serine/threonine kinase 
family (Griner and Kazanietz, 2007), and this kinase is over- 



Figure 3. BRCAl Wild-Type and Mutant 
Expression in Fibroblasts, iPSCs, and Ter- 
atomas 

(A) qPCR for absolute BRCAl mRNA levels in 
wiLd-type (WT) and BRCAl 5382insC fibro- 
blasts and iPSCs. Reprogramming Led to a 
doubling of total BRCAl mRNA in iPSCs. A 
total of eight fibroblast cell Lines and 24 
iPSC lines were analyzed. Each point is an 
average of three replicate PCR reactions. 
The difference in expression of BRCAl be- 
tween the iPSCs and fibroblasts was statis- 
tically significant (Pearson chi squared, p < 
0.001). The western blot for BRCAl (top 
panel) in six BRCAl WT fibroblast (Lanes 
1-3) or iPSC (Lanes 4-6) lines confirmed 
that the BRCAl protein is also induced with 
reprogramming. The HIPl protein was used 
as a loading control and migrates faster 
when derived from iPSC extracts compared 
with fibroblast extracts (bottom panel). 

(B) Allele-specific qPCR was used to quan- 
titate the relative levels of WT (black bar) 
and mutant (gray bar) BRCAl mRNA in 
fibroblasts. The eight fibroblast cell lines 
were analyzed with three replicate reactions 
and averaged. There were no significant 
differences in expression of the WT and 
mutant alleles. 

(C) qPCR, as in (B), was used to quantitate 
the relative levels of WT (black bar) and 
mutant (gray bar) BRCAl mRNA in iPSC 
lines. The 24 iPSC lines were analyzed with 
three replicate reactions and averaged. 
There were no significant differences in 
expression of the WT and mutant alleles. 

(D) qPCR, as in (B), was used to quantitate 
the relative levels of WT (black bar) and 
mutant (gray bar) BRCAl mRNA in tera- 
tomas. RNAs from the 24 teratomas were 
analyzed with three replicate reactions and 
averaged. There were no significant differ- 
ences in expression of the WT and mutant 
alleles. 

See also Figure S4. 

expressed in most gastrointestinal stromal tumors (GISTs) 
(Blay et al., 2004; Ou et al, 2008; Wong and Shelley-Fraser, 
2010) and has recently been found to be elevated in ER- 
negative mouse breast cancers induced by a c-rel transgene 
(Belguise and Sonenshein, 2007). Because PKC-theta has 
been evaluated extensively for expression in GISTs, we 
used the same assay to stain breast cancer samples (n = 
47) and found that PKC-theta was upregulated in more 
than half of these patients' tumors compared with normal 
tissue (see Figure 4E; Table S2). PKC-theta was not ex- 
pressed in normal breast epithelium (n = 39; nine samples 
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Figure 4. Gene-Expression Profiling Identified UpreguLation of PRKCQm BRCAI Mutant iPSCs 

(A) Mi'croarray analysis of gene expression identified PRKCQ, the gene for PKC-theta, as a consistently differentially upregulated gene in 
BRCAI 5382insC iPSC lines (p < 0.001). Twenty-four iPSC lines were analyzed. Eleven were WT and the other 13 were mutant. 

(B) The relative expression levels (blue, zero expression; orange, highest expression) of other members of the PKC family exhibited no 
significant differences between WT and mutant iPSC lines. 

(C) qPCR showed that fibroblast lines express very Low levels of PKC-theta relative to iPSC lines and that PKC-theta is elevated in BRCAI 
mutant iPSCs (p < 0.001). A total of eight fibroblast cell lines and 24 iPSC lines were analyzed. Each point is an average of three replicate 
PCR reactions. The difference in expression of BRCAI between the iPSC lines and fibroblast lines was statistically significant (Pearson chi 
squared, p < 0.001). 

(D) Immunoblot analysis confirmed that PKC-theta protein was also increased in all BRCAI mutant heterozygous iPSC lines. (Top) lanes 
1-6: BRCAI WT iPSC lines (iPS.II.l.l, iPS.II.1.2, iPS.II.1.3, iPS.III.4.1, iPS.III.4.3, and iPS.III.1.1); lanes 7-13: BRCAI mutant iPSC lines 
(iPS.III.7.4, iPS.III.7.1, iPS.III.7.2, iPS.III.7.3, iPS.II.6.1, iPS.II.6.2, and iPS.II.6.3). (Bottom) lanes 1-5: BRCAI WT iPSC lines 
(iPS.III.1.2, iPS.III.1.3, iPS.IV.1.1, iPS.IV.1.2, and iPS.IV.1.3); lanes 6-10: BRCAI mutant iPSC lines (iPS.IV.3.1, iPS.IV.3.2, iPS.IV.3.3, 
iPS.ni.3.1, and iPS.III.3.2). 

(E) PKC-theta immunohistochemistry of representative breast tumors (right-hand panels) that have weak (1), intermediate (2), or strong 
(3) stains. The normal tissue stains weakly (1) for PKC-theta (left). Scale bars represent 50 |j,m. 

(F) qPCR showed that BRCAI mutant breast cancer cell lines (HCC1937 and MDA-MB-436) expressed higher Levels of PKC-theta compared 
with the tested WT breast cancer cell lines. 

See also Table S2. 



were matched to patient tumors), consistent with reports 
of restricted PKC-theta expression in T lymphocytes, plate- 
lets, and skeletal muscle (Isakov and Altman, 2012). Tu- 
mors that expressed high levels of PKC-theta were more 
frequently ER/PR negative than those that did not express 
PKC-theta. And tumors with low PKC-theta expression 
were often ER/PR positive (55%), compared to only 22% 
ER/PR positivity among tumors with high PKC-theta 



expression (Table S2; p < 0.03). No significant correlation 
was found between PKC-theta staining intensity and 
Her2/neu, p53, or tumor stage (Table S2). 

Whether PKC-theta is overexpressed in tumors from the 
subset of breast cancer patients with germline BRCAI mu- 
tations and whether this factor is a predictive marker or 
therapeutic target for some breast cancers is unknown. 
Here, we found that out of five tested breast cancer cell 
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lines, the two breast cell lines known to have deleterious 
C-terminal BRCAl frameshift (HCC1937 has the BRCAl 
5382insC mutation of the family studied here) or trunca- 
tion (MDA-MB-436 has the BRCAl 5396+lG > A mutation) 
mutations had significantly higher PKC-theta mRNA levels 
(Figure 4F). The other three breast cancer cell lines known 
to have wild-type BRCAl (MDA-MB-134-V1, MCF7 and 
T47D) did not express PKC-theta. 

Whole-Genome Sequence Analysis of iPSC Lines 

To determine whether increased expression of PKC-theta 
in BRCAl mutant iPSCs resulted from an accumulation 
of distinct de novo mutations during the reprogramming 
process, we performed WGS on all 24 iPSC lines, their 
parental fibroblasts, and matched lymphoblasts to iden- 
tify de novo mutations. It has been reported that the 
process of reprogramming results in de novo coding muta- 
tions in iPSCs (Gore et al., 2011; Ji et al., 2012; Ruiz et al., 
2013). 

Variants identified by WGS were conservatively analyzed 
to reduce the number of false-positive de novo mutations. 
Because putative de novo variants are often false positives 
with low sequence quality, we saw an enrichment of de 
novo variants with low-quality scores compared with vari- 
ants annotated in the Single Nucleotide Polymorphism 
database (dbSNP; Figure 5A). A stringent quality threshold 
was used to filter de novo variants (see the Experimental 
Procedures). The performance of quality control (QC) 
filtering was assessed using two methods: (1) measuring 
genome-wide genotype concordance between matched 
fibroblast and lymphoblast genomes and (2) Sanger 
sequence validation of randomly selected variants. 
Genome-wide genotype concordance between the fibro- 
blast and lymphoblast genomes after QC filtering was 
>98.45% (data not shown), and Sanger sequencing 
confirmed 8S% of the called variants (Table S3). 

Nonsynonymous sequence variants identified by WGS 
of each iPSC line were analyzed to identify high-confi- 
dence de novo mutations (Table S4). To distinguish de 
novo nonsynonymous coding mutations from those in- 
herited from small populations of parental fibroblasts, we 
excluded iPSC line mutations present in the matched fibro- 
blast and/or lymphoblast sequences. This effectively 
excluded QS.eyo-QQ.So/o of nonsynonymous variants. 
Because the iPSCs in this study were generated from related 
individuals, we further refined our analysis by excluding 
mutations present in iPSCs, fibroblast, and/or lymphoblast 
sequences of other related individuals. This additional 
layer of filtering allows for improved specificity in deter- 
mining whether a variant is a true or false positive call 
(Gore et al., 2011; Ji et al., 2012; Ruiz et al., 2013). It also 
reduces the chance of mistakenly categorizing a false-nega- 
tive variant in the parental fibroblasts as de novo in the 



iPSCs. Because de novo mutations are not likely to appear 
in public databases such as the 1000 Genomes Project and 
the NHLBl Exome Variant Server, we removed all putative 
de novo mutations that appeared in these databases. 
Finally, because de novo mutations are unlikely to occur 
in more than one iPSC line (unless of course they are 
selected for biological reasons), we compared de novo 
mutations between iPSC lines and excluded mutations 
present in more than one cell line. As expected, we 
observed a decrease in the median quality score of variants 
throughout out analysis (Figure S4). Sequence analysis 
identified on average 10.7 nonsynonymous de novo muta- 
tions per iPSC line, ranging from no mutations to 21 
mutations (Table S5), not including the iPSC line 
iPS.IV.3.2, which had a uniquely ~12-fold higher number 
of de novo nonsynonymous mutations. Pathway analysis 
failed to identify any specific enrichment for de novo non- 
synonymous mutations regardless of BRCAl mutation 
status (data not shown), suggesting a random distribution 
of de novo mutations throughout the genome. Consistent 
with this, the frequency of de novo nonsynonymous 
mutations was significantly correlated with the frequency 
of genome-wide de novo mutations for both SNPs and 
insertion/deletions (r = 0.79, p = 0.00003787) and SNPs 
alone (r = 0.82; p = 0.0000164). 

Since heterozygosity for BRCAl mutations has been re- 
ported to result in genomic instability (Konishi et al., 
2011), we expected to observe an increased frequency of 
de novo mutations in the BRCAl mutant iPSCs compared 
with wild-type iPSCs. We compared the number of nonsy- 
nonymous de novo mutations in BRCAl wild-type (n = 9) 
versus mutant (n = 11) iPSC lines to determine whether 
iPSC lines heterozygous for the BRCAl 5382insC mutation 
had an increased frequency of mutations. No statistically 
significant difference was observed in the number of de 
novo nonsynonymous SNPs and insertions/deletions (p = 
0.08053) or SNPs alone (p = 0.06208; Figure 5B; Table S5). 
Likewise, the number of genome-wide de novo SNPs and 
indels or SNPs alone was not significantly different 
between the BRCAl wild-type and mutant iPSC lines (Table 
S6). From these data, we conclude that although de novo 
mutations do occur during the reprogramming of fibro- 
blasts into iPSCs with mRNA transfection technology 
(Warren et al., 2010), BRCAl haploinsufficiency does not 
significantly increase mutation frequency in coding or 
noncoding regions of the iPSC genome. 

One iPSC line (iPS.lV.3.2) displayed all of the normal 
morphological characteristics and differentiation capabil- 
ities assayed in our screen for normal iPSC lines, but 
acquired -^12 times as many de novo mutations as all other 
iPSC lines (Tables S5 and S6). Although the cause of this 
genomic instability is not known, a de novo nonsense 
mutation in PALB2 (Q1004X) was identified by WGS and 
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Figure 5. Whole-Genome Sequence Analysis of iPSC Lines 

(A) Distribution of quality scores for genome-wide SNPs annotated in the dbSNP database (black) compared with de novo SNPs (yellow). 
The y axis "proportion" is the frequency of all identified variants (numerical distribution). Each line represents a single genome's 

(legend continued on next page) 
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confirmed by Sanger sequencing. PALB2 binds to the 
coiled-coil domain of BRCAl and collaborates with 
BRCAl in homologous recombination DNA repair 
(Sy et al., 2009). Additionally germline PALB2 mutations 
have been implicated in familial pancreatic cancer syn- 
dromes (Jones et al., 2009). It remains possible that muta- 
tions in both BRCAl and PALB2 lead to increased genome 
instability. In addition to an increased frequency of de 
novo mutations, iPS.IV.3.2 had de novo structural abnor- 
malities. WGS coverage analysis suggested copy-number 
variation of chromosome 20, including loss and gain of 
copy number of the entire short (p-arm) and long (q-arm) 
arms, respectively LOH analysis confirmed de novo LOH 
for chromosome 20p (Figure 5C). Consistent with this, 
expression analysis showed expected changes in gene 
expression based on the copy-number variations of genes 
on the p-arm of chromosome 20 compared with other 
iPSC lines from patient IV.3 (Figure 5D). 

DISCUSSION 

We report the generation of a series of iPSC lines from indi- 
viduals with the BRCAl 5382insC mutation, which is a 
well-known predictor of high risk for breast and ovarian 
cancer (Futreal et al., 1994). Generation of iPSCs from indi- 
viduals with BRCAl mutations provides an opportunity to 
examine why cancer develops more frequently in these pa- 
tients. In this study, it was easier to reprogram the mutant 
fibroblasts than the fibroblasts from relatives without the 
BRCAl mutation. Determining whether this difference is 
due to the mutation itself or simply to line-to-line varia- 
tions due to other characteristics (e.g., growth characteris- 
tics and handling) will require the generation of iPSCs 
from additional families and/or heterozygous mice with 
BRCAl mutations. Our data contrast with those obtained 
by Gonzalez et al. (2013) from mouse embryonic fibro- 
blasts (MEFs), in which homozygous mutations of Brcal 
led to reprogramming defects. Those authors did not 



include MEFs from heterozygous Brcal mutant embryos, 
and the human BRCAl mutations we used are different 
from the mouse Brcal mutations used in their study. 
Finally, the parental fibroblasts in the mouse experiments 
were embryonic, whereas our fibroblasts were adult. The re- 
programming experiments in mice with the BRCAl 
5382insC mutation reported here, and TALEN-mediated 
targeting of the remaining human wild-type BRCAl allele 
in the heterozygous iPSC or fibroblast lines generated in 
this study will begin to address these differences. 

We observed that mutant and wild-type BRCAl alleles 
were expressed at equivalent levels in fibroblasts and their 
derivative iPSCs and teratomas. There was also a specific 
upregulation of PKC-theta in the BRCAl 5382insC iPSCs 
compared with the wild-type iPSCs. A functional relation- 
ship between BRCAl and PKC-theta levels may involve 
the known interaction of BRCAl with the HDAC2 complex 
(Chang et al., 201 1) or the role of BRCAl as a transcription 
factor (Haile and Parvin, 1999; Horwitz et al., 2006; Schle- 
gel et al., 2000). 

A related finding in this report is that PKC-theta is over- 
expressed in some breast tumors. PKC-theta has been 
found to be a reliable marker of GISTs (Blay et al., 2004; 
Motegi et al., 2005; Ou et al., 2008), and prior reports 
have suggested that PKC activity, although not specifically 
PKC-theta activity, is upregulated in a subset of breast can- 
cers (Gordge et al., 1996; O'Brian et al., 1989). We found 
PKC-theta to be upregulated in a fraction of breast tumors 
that are frequently hormone-receptor negative, a charac- 
teristic of BRCAl mutant tumors. Because of the tumori- 
genic potential of PKC-theta, these results suggest that 
lower levels of PKC-theta in the presence of wild-type 
BRCAl may be a mechanism for the tumor-suppressor 
function of BRCAl. In fact, Belguise et al. determined 
that PKC-theta is elevated in ER-negative human breast 
cancer cell lines (Belguise et al., 2012) and NF-kappaB- 
induced mouse breast tumors (Belguise and Sonenshein, 
2007). They showed that PKC-theta inhibits ER expression. 
The MDA-436 BRCAl mutant cell line, which expresses 



distribution (iPSC, fibroblast, or lymphoblast Line). De novo variants are generally called with low quality. A quality score threshold of 100 
(blue line) was used to exclude potential false-positive sequencing errors called at low quality. 

(B) Frequency of de novo nonsynonymous mutations in iPSC lines with WT and mutant BRCAl. No significant difference in the frequency of 
de novo nonsynonymous mutations was observed between WT (n = 9) and mutant (n = 11) BRCAl iPSC lines. The iPS.IV.3.2 outlier was 
excluded from this analysis. "All" mutations include SNPs and indels (left-hand side) and"SNPs only" (right-hand side) are displayed as bars 
that represent the numbers of de novo variants identified per cell line. We calculated p values to detect differences in mutation frequency 
between mutant and WT using Student's unpaired ttest. No significant differences were identified (p <0.08 and p < 0.06 for "All" and "SNPs 
only," respectively). 

(C) De novo LOH of chromosome 20p in iPS.IV.3.2 (upper panel) as determined from the average heterozygosity of SNPs using a 20 KB 
segmented window. No LOH was detected in the matched fibroblast sample (lower panel). 

(D) Copy-number variation in chromosome 20p of iPS.IV.3.2 was confirmed by transcriptional profiling of genes on this chromosome. 
Heatmaps show the fold change in expression in iPS.IV.3.1, iPS.IV.3.2, and iPS.IV.3.3 versus the average expression level in the other 21 
iPSC lines. Red, <0. 5-fold; green, >2. 0-fold. 

See also Figure S5, Table S3, and Tables S4-S6. 
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PKC-theta at levels higher than any of the other ER-nega- 
tive breast cancer lines tested in their report, was the only 
BRCAI mutant cell line they analyzed (Belguisc et al., 
2012). However, all of the ER-negative lines that were 
found to have high PKC-theta have been reported to be 
haploinsufficientforBJ^C^J (Elstrodt et al., 2006). Upregu- 
lation of PKC-theta in our BRCAI mutant iPSCs may be 
direct or indirect. If BRCAI is directly upstream of PKC- 
theta, this would explain why BRCAI mutant breast can- 
cers are frequently ER negative (Hosey et al., 2007). 

Although our work with BRCAI mutant iPSCs led to our 
investigation of PKC-theta upregulation in human breast 
tumors, future studies will be necessary to validate whether 
PKC-theta can be used as a biomarker or therapeutic target 
in BRCAI mutant tumors. Because the frequency of BRCAI 
mutations in patients with breast cancer is only 5%-10%, 
but we found PKC-theta to be elevated in greater than 
half of the breast tumors in this study, the PKC-theta 
upregulation must occur in some tumors without BRCAI 
mutations. PKC-theta elevation does support the consider- 
ation of this protein as a new drug target, with potential 
applications in patients with BRCAI mutations (Figure 4E; 
Table SI). 

Through our use of WGS, we increased the depth of our 
characterization of iPSCs to include not just cellular but 
also genetic changes that occur during the reprogramming 
process in BRCAI mutant cells. It has been reported that re- 
programming of differentiated cells to a pluripotent state 
results in an estimated 6 protein-coding point mutations 
per exome (Gore et al., 2011). By comparing the WGS of 
iPSCs generated in this study with not only their fibroblast 
progenitors but also the fibroblasts of related individuals, 
we increased the ability to remove mutations that are not 
de novo but instead are inherited from the fibroblast pro- 
genitors. To our surprise, the WGS approach identified 
only slightly higher numbers of de novo nonsynonymous 
mutations than Gore et al., although the small difference 
can be attributed to the difficulty of calling high-confi- 
dence insertion/deletion mutations using hybridization 
capture/sequencing methods (Gore et al., 2011). If we 
consider only SNP mutations, our results (average 6.7 non- 
synonymous de novo SNPs per iPSC line) are consistent 
with previous reports (Gore et al., 2011). To our further 
surprise, through this analysis we did not find a statistically 
significant difference in the number of de novo nonsynon- 
ymous mutations in BRCAI wild-type versus mutant iPSC 
lines. This contrasts with the prior report that the human 
BRCAI heterozygous state is genetically unstable (Konishi 
et al., 2011). It is intriguing to speculate about why our 
result differs from theirs. Konishi et al. analyzed genome 
stability in breast epithelial cells, which are known to be 
susceptible to BRCAI mutant tumorigenesis. Because our 
genome analyses were conducted in stem cells, it will be 



important to analyze the genomes of their derivative 
epithelial cells. One exception to the surprising lack of 
genome instability in BRCAI mutant iPSC lines was the 
mutant iPSC line iPS.lV.3.2, which harbored >10-fold 
more de novo mutations than the other iPSC lines. 
Although this resulted in significant differences in the 
gene-expression profiles, we could not readily identify a 
functional consequence of this genomic instability using 
standard iPSC characterization techniques. This indicates 
that WGS of iPSC lines prior to biological study is 
warranted. 

In conclusion, we have generated and characterized 
patient-derived iPSC lines with a BRCAI mutation and 
identified a correlation between BRCAI mutations and 
PKC-theta expression. Similar analyses of other mutant 
iPSC lines may lead to the discovery of other molecular 
consequences that result from specific mutations associ- 
ated with BRCAI. These data should encourage others to 
generate iPSC lines from families with mutations that 
confer a high cancer risk. 

EXPERIMENTAL PROCEDURES 

Generation and Characterization of Human iPSCs 

Primary fibroblasts were established from skin biopsies taken from 
wild-type and BRCAI 5382insC heterozygous mutant family mem- 
bers shown in Figure lA and Table 1 after informed consent and 
with institutional review board approval. Skin-punch (forearm) 
biopsies were cultured in six-well plates at 37°C in 5% CO2 in 
1 ml Dulbecco's modified Eagle's medium (DMEM) with 10% fetal 
bovine serum, penicillin, and streptomycin until fibroblasts 
migrated out from the skin explants. The iPSCs were generated 
by nonintegrative mRNA-based reprogramming using a cocktail 
of mRNAs (and in some cases miRNAs) for 0CT4, S0X2, KLF4, 
LIN28, and c-MYC as previously described (Anokye-Danso et al., 
2011; Miyoshi et al., 2011; Warren et al., 2010). All 24 cell lines 
are available for research purposes. 

In Vitro and In Vivo Differentiation into Teratomas 

iPSCs were subjected to directed differentiation in vitro as previ- 
ously described (Chambers et al., 2009; D'Amour et al., 2005, 
2006). For the teratoma formation assays, iPSCs (approximately 
1-2 X 10^ cells) were injected into the kidney capsules of immuno- 
deficient NOD/SCID mice. Teratomas were harvested 10-16 weeks 
later. A small portion from each teratoma was extracted for RNA 
analysis and the rest was fixed in PBS containing 4% paraformalde- 
hyde for histology. 

mRNA Expression Profiling in Fibroblasts and iPSCs 

For gene-expression analysis in primary fibroblasts and 
iPSC lines, labeled complementary RNA was hybridized to an 
HT-12v4 BeadChip (Illumina). Data analysis was performed 
via fold-change analysis of relative signal intensities followed 
by Student's paired t test. A p value < 0.05 was considered 
significant. 
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RNA Isolation and qPCR Analysis of Gene Expression 

RNA was purified with Trizol (Invitrogen) and treated with 
DNAase (QIAGEN) to remove genomic DNA contamination. 
One microgram of total RNA was used for the reverse transcriptase 
reaction using Script III (Invitrogen) and random hexamer primers 
according to manufacturer's instructions. qPCR was performed us- 
ing the StepOnePIus Real-Time PGR system (Applied Biosystems) 
with the Taqman gene-expression master mix. Mismatched 
primer sequences for BRCAl gene analysis are described in Fig- 
ure SI. Standard Taqman primers for PRKCQ were purchased 
from Invitrogen. 

Immunohistochemistry 

Tissue microarrays (IMH-364 and IMH-371) of breast cancers and 
normal adjacent tissues were purchased from Imgenex and sub- 
jected to staining with anti-PKC-theta antibody (mAb 27; BD 
Bioscience) as previously described (Motegi et al., 2005). Scoring 
was performed by two scorers (D.R. and T.S.R.) using intensity 
scores of 1-3 (1, weak staining; 2, intermediate staining; 3, strong 
staining). 

Western Blot Analysis 

The iPSGs were homogenized in Ix RIPA buffer (Gell Signaling 
Technology) containing 1 mM phenylmethylsulfonyl fluoride, 
centrifuged at 10,000 rpm for 15 min at 4°C, and protein quanti- 
tated (DC assay; BioRad). Eighty micrograms of total protein was 
resolved on SDS-PAGE, transferred to nitrocellulose, and blotted 
for PKC-theta (mAb 27; BD Biosciences), a-tubulin, anti-HIPl 
4B10 monoclonal, and BRGAl monoclonal (Abeam). 

WGS and De Novo Mutation Identification 

WGS using the combinatorial probe-anchor ligation method 
was performed by Complete Genomics as previously described 
(Drmanac et al., 2010). WGS of all iPSC lines, fibroblasts, and 
lymphoblasts generated sufficient coverage for variant calling 
and copy-number analysis, both performed by Complete Geno- 
mics. QG was performed by excluding iPS variants if (1) the 
variant position was "no-called" in two or more samples (to re- 
move genomic regions that would be difficult to assemble using 
short-read technology) or (2) the variant quality score was <100. 
QG was confirmed by measuring concordance between matched 
fibroblast and lymphoblast samples (>98.45) and direct Sanger 
sequencing of randomly selected mutations (85% confirmed). 
Nonsynonymous iPSG variants passing the QG filters were 
analyzed to identify de novo mutations by excluding (1) variants 
in the fibroblast or lymphoblast of sample-matched or related 
family members, (2) variants annotated in the Exome Variant 
Server containing ~6,500 exomes or the 1000 Genomes database, 
or (3) variants present in any of 54 samples sequenced by Com- 
plete Genomics or whole-exome sequenced using the Illumina 
platform (Table S3). Finally, similarly to previous reports (Gore 
et al., 2011), only variants present in single iPSG lines were 
considered de novo. 

LOH analysis of SNPs passing the QG filters was performed on 
iPS.IV.3.2 and the matched fibroblast (Fb.IV.3). Average heterozy- 
gosity was calculated using a 20 KB segmented window and plotted 
along the entire chromosome 20. 



Confirmation of De Novo Mutations by Sanger 
Sequencing 

For each iPSG line, a subset of the putative de novo mutations was 
selected for confirmation by Sanger sequencing. DNA was ex- 
tracted from iPSC and fibroblast cell lines by standard procedures. 
Primers were designed using Primer Express 3.0 software to amplify 
the genomic region in which each mutation was located. The 
cycling conditions were as follows: 4 min at 94°G, 40 cycles of 
1 min at 94° G, 30 s at the appropriate annealing temperature for 
each primer (49-54°G), 1 min at 72°C, and 72°G for 5 min. PGR 
products were resolved on a l%-2% agarose gel. PGR products 
were purified using the QiaQuick Gel Extraction Kit (QIAGEN) 
and sequenced. 

SUPPLEMENTAL INFORMATION 

Supplemental Information includes five figures and six tables and 
can be found with this article online at http://dx.doi.org/10.1016/ 
j.stemcr2013.08.004. 
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